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Abstract 

A variety of viruses tightly pack their genetic material into protein capsids that are barely large enough to enclose 
the genome. In particular, in bacteriophages, forces as high as 60 pN are encountered during packaging and ejection, 
produced by DNA bending elasticity and self-interactions. The high forces are believed to be important for the 
ejection process, though the extent of their involvement is not yet clear. As a result, there is a need for quantitative 
models and experiments that reveal the nature of the forces relevant to DNA ejection. Here we report measurements 
of the ejection forces for two different mutants of bacteriophage A, Ab221cI26 and AcI60, which differ in genome 
length by ~30%. As expected for a force-driven ejection mechanism, the osmotic pressure at which DNA release is 
completely inhibited varies with the genome length: we find inhibition pressures of 15 atm and 25 atm, respectively, 
values that are in agreement with our theoretical calculations. 

1 Introduction 

Over the past thirty years, a series of experiments and theoretical work have produced many insights about the impor- 
tance of internal force s in the bacte riophage li fe cycl e: Early measurements on A capsids showed that they contained 
tightly packed DNA fearnshaw and Harrisonl Il977l) . and subsequent experiments established that DNA packed at 
these densities exerts a pressure of tens of atmospheres that is dependent on the density and salt conditions jRau et ali 
Il984l) . Any effect of the A genome length on its life cycle (independent of any particular genes) suggests that in- 
ternal forces are important, and there are several such effec ts known: t here are upper and lower bounds on the 
amount of DNA that can be packaged into a A c apsid ( iFeiss et al I ll977l) : mutants of A with long genomes fail to 
grow without magnesium ions jArber et q/.lll983l) : and mutants with short genomes fail to infect pel~ cells (Katsura, 
Il983l) . While magnesium ions reduce the for ces between DN A, stabilizing the phage particles, DNA-condensing 
ions su ch as putrescine pre vent DNA ejection jKatsuraL Il983l) . and osmotic stress can stabilize the genome within 
phages dSerwer et a/.lll983l) . The evidence seems to indicate that internal forces have an important role in the function 
of A: phages with low forces can be incapable of ejecting their DNA forcefully enough to penetrate the cell, while 
phages with high forces are unable to package their genome or are unst able when fully packaged. A variety of theo- 
retical models of DNA packaging in bacteriophage h ave been proposed (Riemer and Bloomfield, 1978; Black, 1989; 
ISerwedJ^SStlTzhT et al. l l2Q03UPurohit et al\XMS% b ut only recently have experiments begun to quantify the forces 
required to tightly pack DNA into capsids ( Smit hef all 120011) and the forces driving DNA ejection dEvilevitch et all 
20031 120051) . 
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Figure 1: Comparison of Ab221 (lane 2) and AcI60 (lane 3) genomes by field-inversion gel electrophoresis. The gel 
was run with 100/60 V switching for 19h in 0.5x TBE buffer, 1% agarose. Lanes 1 and 4 are A-mix ladders, with 
known lengths as shown. Bands contain 0.5-1 ng of DNA, stained with SYBR Gold. The gel shows a single band in 
lane 2 at 37.9 ± 0.3 kbp and a single band in lane 3 at 48.4 ± 0.3 kbp; both results are consistent with the expected 
values of 37.7 kbp and 48.5 kbp. 

The aim of this paper is to study the effect of genome length on the ejection force of A DNA, by comparing AcI60, 
a simple mutant of the wild ty pe with a 48.5 kbp genome, to Ab221cI26 (Ab221), which has a much shorter genome of 
37.7 kbp (Bell ett et aZ.LIl971l> . The reason that measurements with different genome lengths are especially interesting 
is that simple models of the forces that arise during packaging depend in a precise way upon the genome length. To 
measure the force, a method reported earlier jEvilevitch et all 120031) was used: osmotic stress was applied to the 
outside of the capsids during ejection, halting the ejection at the point where the internal and external forces balance. 
Using this method, we show that phages with shorter genomes have lower forces, phages with longer genomes eject 
their DNA with higher forces, and that straightforward theoretical models are sufficient to predict these effects. 

2 Materials and Methods 

Phages Ab221cI26 (Ab221) and AcI60 were extracted from single plaques, grown in 3 L cultures of E. coli c600 cells, 
and purified by PEG precipitation, differential sedimentation, and equilibrium CsCl gradients, resulting in ~10 13 
infectious particles. After purification, phages were dialyzed twice against a 500-fold greater volume of TM buffer 
(50 mM Tris, 10 mM MgS0 4 , pH 7.4). 

To check the genome lengths of the phages used in this experiment, we extracted the DNA with phenol and 
chloroform from approximately 5 x 10 9 phages of each type into 500 /iL of 0.5x TBE buffer. A quantity of 1 /xL 
(10 7 genomes, or 0.5 ng) was removed from each extraction, mixed with loading dye, heated briefly to 65°C to 
separate cohesive ends, and pipetted into a 0.5x TBE, 1% agarose gel. A 10 ng quantit y of a sta ndard ladder (A-mix, 
Fermentas) was included for size comparison. We ran the gel using the method of ( Birr en et al ., 1990), with a simple 
electrophoresis box (Owl Separation Systems B1A) and a homebuilt voltage inverter, pulsed at 100 V forward, 60 V 
backward, for 19 h. The gel was stained with SYBR Gold (Molecular Probes) and photographed with an Alphaimager 
HP (Alpha Innotech). Dust was manually erased from the resulting image, which is shown in Fig. \l\ Results were 
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Figure 2: UV absorbance of the supernatant containing fragments of the DNA ejected from Ab221 (data set "Ab221 
set 1", left) and AcI60 (data set "AcI60 set 1", right), aligned at 300 nm. Absorbance curves are shown for various 
concentrations of PEG, corresponding to different osmotic pressures, and for tests done in the absence of LamB. A 
second data set taken for each phage is not shown. 



consistent with the expected 37.7 kbp genome for Ab221 and 48.5 kbp genome for AcI60. 

The A receptor LamB (maltoporin), required to trigger ejection, was extracted from the membranes of E. coli 
popl54 cells. These cells express a lamB gene from S. sonnei known to be compatible with a variety of A strains, 
allowing ejection in the absence of chloroform (Roa and Sc andellal Il97<4 iGraffef a/.L |2002). LamB was affinity- 
purified in amylose resin and dialyzed twice against TM buffer containing 1% n-octyl-oligo-oxyethylene (oPOE; 
Alexis Biochemicals #500-002-L005.) 

Our method for measuring ejection forces is substantially the same as that described earlier jEvilevitch et all 

l2003h . with minor refinements that have improved precision. We paid particular attention to the difficulty of pipetting 
the viscous PEG solutions, trying to minimize systematic and statistical errors that occur when the solution adheres 
to the pipette tips. A solution of 50% (w/w) polyethylene glycol (PEG) 8000 (Fluka PE G Ultra) was prep ared in 
TM buffer with 0.5% oPOE, and its density was measured at 1.09 g/ml (see also ( Gonzalez-Tell o et al l ll994l) ). This 
solution was used to prepare solutions of PEG / TM 0.5-1% oPOE at various specified %(w/v) values (see Fig. [2) on 
an analytical balance. The mass measurements allowed us to set the quantity of PEG in each sample within 0.2 mg, 
which corresponds to an error in concentration of approximately 0.1% (w/v). Phage solution was added to a final 
concentration of ~10 1:l /ml, and DNase I was added at 10 /ig/ml. The sample tubes were turned slowly for several 
minutes to mix the viscous PEG solutions. Purified LamB was added with a wide-mouth tip and the resulting 200 fxL 
solution was mixed quickly by pipetting. A wide-mouth tip is an inaccurate pipetting device, but it is necessary for 
quickly mixing viscous solutions. To minimize the effect of inaccurate pipetting on the measurements, we used a 
concentration of LamB that was sufficient for maximal ejection, estimated at ^5 /ig/ml by UV absorbance. After 
the addition of LamB, the samples were incubated for 1 h at 37°C, which was sufficient for the reaction to reach its 
endpoint — complete digestion of the ejected genome fraction by DNase I. Finally, the capsids were separated from the 
ejected DNA fragments by a centrifugation for 20 h at 18,000 x g. 

After centrifugation, 1 20 [iL of supernatant from each tube was removed to a U V-transparent plastic cuvette (Ocean 
Optics UVettes) and DNA concentrations were measured with a UV-visible spectrophotometer (LKB Biochrom Ultro- 
spec II). The absorbance curves were aligned at 300 nm (immediately after the DNA absorbance peak) to compensate 
for absorbance not due to DNA. The resulting curves are shown in Fig. [2] The absorbance values at 260 nm, A260> are 
linearly related to the amount of DNA ejected from the phage capsids. In contrast with earlier experiments, there was 
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no measurable background DNA absorbance due to ruptured phage capsids: samples at high PEG concentrations or 
without LamB had similar values of ^260 to samples prepared without phages. This is probably because phages were 
used within one month of dialysis; in contrast, samples of Ab221 measured after five months of storage in TM buffer 
at 4°C had a background A 2 eo ~ 0.1. 

When there is no PEG present, ejection reaches completion (100%), and when no phages are added there is no 
ejection (0%). Intermediate values were found with a linear interpolation: 

ejected fraction = 100% • (A 260 with PEG - ^260,no phages)/ ( A 260,no PEG ~ ^260,no phages) ■ (*) 

An alternative procedure is to use no LamB as a calibration for 0% (Evi levitch et a/l E005). 

The weight measurements set ttipeg, the mass of PEG, and mtot, the total mass, for each sample. The % (w/w) 
weight-weight fraction was computed a s w = mpreo/ mtnt.. The osmotic pressure at each PEG concentration was then 
determined with the empirical formula (Michel, 1983) 

n(atm) = -1.29G 2 T+ 140G 2 + 4G, (2) 

where T is the temperature (°C) and G = w/(100 — w). Note that the osmotic pressure is an increasing function of 
the PEG concentration and a decreasing function of temperature. For this experiment, T = 37. 

Two complete sets of samples were prepared for each of AcI60 and Ab221. Each sample has a statistical error due 
to weight measurements, pipetting, and spectrophotometry. To minimize systematic effects on the ejected fraction, the 
two "no phage" and "no PEG" tubes were averaged for each phage. Statistical errors were propagated to yield x and 
y error bars (see Fig.|3|) 



3 Theoretical Model of DNA Packaging 

Our t heoretical model is based on earlier work dRiemer and BloomfieldLll978l:lTzril et a/.LE003llPurohit et q/.ll2003l 
120051) . which describes the packaging energy as a function of the length of DNA in the capsid. We model the A 
capsid as a sphere and its genome as a long semiflexible rod. We assume that the rod is wound into a cylindrically 
symmetric spool with local hexagonal packing. The total energy of the packaged DNA can 

then be approximated by a sum of inter-axial repulsion energy and the bending energy of the rod: 

(■-Rout jyfjA 

E = ^interaction + £Wid = V3F L(c 2 + cd s ) exp(-d s /c) + 7rfc B T£ / dr , (3) 

where Fo and c are empirically determined constants describing the interaction between neighboring DNA double- 
helices, £ is the persistence length of DNA, L is the length of the DNA within the capsid, d s is the inter-axial spacing, 
i? ou t and i?; n are the radius of the capsid and the inner radius of the DNA spool, respectively, and N(r) is the number 
of loops of DNA at a distance r from t he spool axis. For the persistence length £ we use 50 nm, though its value 
in Mg 2+ buffer may be ~10% smaller (Hagerman, 198 oj. The spacin g between sequential bases of DNA varies, 
depending on the base types, from 0.33 to 0.34 nm dOlson et all fl998). To compute L we disregard this variation 
and use 0.34 nm tim es the number of base pairs within the capsid. The inter-axial forces in buffers containing Mg 2+ 
have been measured dRauefadfl984l) . Since the values measured for 5 mM and 25 mM Mg 2+ were not significantly 
different, we assume that t he forces at 10 m M (used in our experiments) will be identical. A least-squares fit to the 
5 mM and 25 mM data in dRau et ad 1 1984 gives F n = 12,0 00 pN/nm 2 and c = 0.30 nm. The radius of the phage 
capsid R ou t is around 29 nm dEarnshaw and Harrison! 1 19771) . Once we know d s , i?; n , and N(r), we can use Eq. [3] 
to calculate the internal force on the phage genome as a function of genome length inside the capsid, providing an 
interpretation of the experimental results. 

We calculate the remaining variables as a function of L according to the following recipe, which involves only 
simple geometrical considerations and elementary calculus. The number of loops N(r) in Eq. |3]is given by z(r) /d s , 
where z(r) = yj i? 2 ut — r 2 is the height of the capsid at distance r from the central axis of the DNA spool. The actual 
volume V(Ri n , R ou t) occupied by the DNA spool can be related to the genom e length L and the inter-axial spacing 
d s to get an expression for R m in terms of d s , R ou t and L dPurohit et artl2003l) . This expression for i?; n is substituted 
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into Eq. [3] which can then be minimized with respect to d s to give the equilibrium inter-axial spacing as a function 
of the genome length L inside the capsid. From d s , Rin, and N(r), Eq.[3]now gives us the total packing energy as a 
function of genome length inside the capsi d. The interna l force F(L) acting on the genome is obtained by taking the 
derivative of Eq.[3]with respect to L dPurohit et al 

The preceding construct is a parameter-free model that predicts the ejection force from a A capsid. Experimental 
uncertainties in the quantities quoted above should lead to errors of 10-50% in the magnitude of the force predicted, 
but the shape and relative positions of the curves for different genome lengths should not be strongly affected by these 
errors (these tests of the parameters are not shown). 

DNA ejection in our experiment or in vivo is halted at the point where the internal force balances the osmotic 
force. We have described above how to obtain the internal force (and hence, the external osmotic force beca use of the 
equilibrium) acting on the genome. But the experimental variable is an osmotic pressure dTzlilef aZlEool) . Thus we 
need to translate this force into a pressure. The force F(L) is given approximately by IT • 7ri?p NA , where IT is the 
osmotic pressure and -Rdna is the effective radius of the DN A. We take -Rtina as 1.0 nm (bare DNA) plus 0.2 nm, 
half the PEG monomer length found experimentally dAbbot et al.U 19921 iMarshi 120041) . This is an exact formula for 
the osmotic force on a very large area. However, at the PEG concentrations u sed in our e xperiment, the diameter 
of DNA is comparable to the correlation length (mesh size) of PEG, ~l-3 nm ddeGennesifl97i p p.78-80) so the 
formul a for F(L) is only approximately valid. A scaling expression for the correction could be used dde VrieslEoOll 
ICastelnovo et all l2003tlEvilevitch et all 120041) . but since this result is good only up to a multiplicative constant, its 
importance is unclear. Hence, to maintain clarity in our analysis we do not use the correction. On the other hand, 
including this effect could result in a better fit between theory and experiment. 



4 Results and Discussion 

Fig-0shows our experimental results. Partial ejection in both Ab221 (37.7 kbp) and AcI60 (48.5 kbp) corresponds 
closely to that measured earlier in EMBL3 (41.5 kbp): the ejected fraction decr eases quickly to 50% a t several atmo- 
spheres of osmotic pressure, then descends more slowly to 0%. (As shown bv lEvilevitch et al\ d2005h . this "ejected 
fraction" corresponds to all phages ejecting the same percentage of their DNA, rather than a percentage of the phages 
ejecting all of their DNA.) Most importantly, at every pressure where a measurable amount of DNA is ejected, Ab221 
ejects less DNA than AcI60. Conversely, at any point in the ejection process, the osmotic force on the DNA is always 
lower in Ab221 than in AcI60. For example, the highest forces occur at the point where osmotic pressure completely 
inhibits the ejection of DNA. This inhibition pressure is in the range 10-15 atm for Ab221 and 20-25 atm for AcI60 
— a ^100% increase in pressure results from a 30% increase in genome length. 

An interesting dip (or bump) in the data is apparent between 1 .5 and 3.4 atm. We observed this dip in all data sets 
collected for this experiment, but its cause is unknown. 

Fig-0also shows the predictions of our parameter-free theory. Except for the dip, the theory predicts the data quite 
well. Both the absolute magnitude and the general shape of the curves are predicted correctly. This is a remarkable 
result because no fitting was done to match the theory to the data. A straightforward inverse-spool model of the DNA 
arrangement, taking into account only the measured bending elasticity and interaxial forces, correctly predicts the 
dependence of ejection force on genome length. 

Our measurements can be compared with previous results on EMBL3 dEvilevitch et all l2003h . for which the 
genome length is 41.5 kbp. To do so, however, requires that the weight fractions be corrected in calculating the 
osmotic pressures, which we have done here for the rest of the data. In this earlier work, the osmotic pressures of the 
solutions were calculated from the %(w/v) concentrations, without taking into account the density of the PEG solution. 
Here, we convert it to %(w/w) taking into account the density of the PEG solution. The the density differences can be 
corrected using the relation between weight fraction and density of PEG solutions dGonzalez-Tello et aZ.lll994T) . After 
this correction, we find that the maximum pressure, for example, is 15 atm instead of 19.6 atm. The correction has a 
smaller effect at lower PEG concentrations; the difference at 10% (w/v) is only 0.2 atm. Fig.[3]shows a comparison 
between EMBL3 and the mutants used in the present study. Within the experimental uncertainty, the ejection curve 
for the intermediate length genome lies between those of the 48.5 kbp and 37.7 kbp phages reported here. 

An alternative view of the ejection data is shown in Fig.|4] From the ejected fraction p plotted in Fig. [3] we can 
compute the amount of DNA remaining in the capsid as N ■ (1 — p), where N is the number of base pairs in the 
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Figure 3: Measured ejection percentages for four sets of data on AcI60 and Ab221 (points) compared to theoretical 
predictions (curves) for 37.7 kbp and 48.5 kbp genomes. Error bars showing one standard deviation were computed 
from weight measurement, pipetting, and spectrophotometer errors. At all pressures below 15 atm, AcI60 ejects more 
DNA than Ab221. There is an interesting "dip" in the experimental data for both phages, between 1.5 and 3.5 atm, but 
otherwise the precise trends predicted by theory are seen in the data. The theoretical pressures are calculated with the 
approximate relation II (atm) = Ff ttR^ j a , where i?r>NA is the e ffective radius of DNA, assumed to be 1 .2 nm. Also 
plotted: data reported earlier on EMBL3 dEvilevitch et ali 120031) . with corrections made for a systematic error in the 
earlier analysis. There is a larger scattering error in the EMBL3 data than in the AcI60 and Ab221 data; taking this into 
consideration, the data sets for all three phages are completely consistent with the expectation that phages with longer 
genomes eject more DNA at all pressures. 
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Figure 4: Measured ejection data, plotted as the amount of DNA remaining in the capsid after ejection. The remaining 
DNA is calculated as N ■ (1 — p), where N is the number of base pairs in the genome and p is the ejected fraction 
plotted in Fig. [3] Dotted lines at 48.5, 41.5, and 37.7 kbp correspond to the full genomes of AcI60, Ab221, and EMBL3, 
respectively. As expected, the data points fall on a single curve below 10 atm, then diverge as the maximum pressure 
is reached for each phage. 



genome. Since ejected DNA is digested by DNase I, the amount remaining in the capsid should depend only on the 
external osmotic pressure, and it should be independent of the genome size of the phage. Indeed, Fig.|4]shows that all 
of the data below 15 atm falls roughly on a single curve. Above 25 atm, the pressure is sufficient to hold the entire 
genome within the capsid, so the amount of encapsidated DNA is just equal to the genome length of the phage. 

For AcI60, which is almost identical to the wild-type A, our model and measurements indicate that DNA is ejected 
with a force of around 10 pN that drops steadily as the genome enters the host cell. What advantage does A gain from 
having such a high ejecti on force? For the phage DNA to enter the E. coli cell, it must overcome an internal osmotic 
pressure of 3 atm (Neidhardt, 1996). Fig.[5]shows that AcI60 can eject roughly 60% of its genome, 30 kbp, against this 
osmotic pressure. However, the rest of the genome must be actively transported into the cell. Reversible diffusion, in 
particular, is incapable of transpor ting the rest of the genome int o cell, since the energy barrier that must be overcome 
is several times 10 3 kBT ( Smith et al. , 2001 ; Purohit et al. , 2005). A two-step process is implied: first a quick pressure- 
driven injection of half of the DNA, then a slower protein-driven importation of the remainder (as, f or example, by 
RNA polymera se in phage T7). This two-st ep process has been observed in vivo for the phages T5 (Letellier et al., 
|2004 and 029 dGonzalez-Huici et aZ.lEobH) . though similar experiments have not been done on phage A. For Ab221, 
Fig.[3]shows that only 40% of the genome, 15 kbp, is injected in the first step. Since phages with shorter genomes 
than AcI60 are not infectious, we speculate that 15 kbp is close to the minimum amount of DNA required to initiate 
the second (non-pressure-driven) stage of importation. 

In the case of bacterio phage T7, however, the entire genome is actively transported at a slow, constant speed ( M olineu: 
l200lUKemt) ^Zll2004 . even though calculations indicate its internal pressure to be similar to that in A (Pur ohit et all 
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120051) . Apparently, something unique to T7 slows down the process (perhaps the observed expulsion of internal pro- 
teins from the T7 capsid?), whereas most ds DNA bacteriophages experience a first stage of fast, p ressure-driven 
ejection. Since ejection from A is much faster INovick and Baldeschwieler, 1988; Ke mp et a/jj2 004) than fr om T7 , 
it s eems that A is in t his mo st common category. (As another example, recent experiments by d e Frutos et al.\ |2005) 
and [Mangenot et at show that DNA ejection from bacteriophage T5 proceeds very rapidly between discrete 

stopping points along the genome. Furthermore, it has been shown in unpublished measurements that T5 ejection can 
be inhibited by external osmotic pressure.) In vivo ejection experiments with A and in vitro experiments with T7 need 
to be done to explore the differences between these two phages. 
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